Antiferromagnetic Interactions and the Superconducting Gap Function 
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Spin-fluctuation- mediated superconductivity is conventionally associated with d x 2_ y 2 pairing. We 
show that a generalized model of antiferromagnetic spin fluctuations in three dimensions may also 
yield a state with formal "s-wave" (Ai B ) symmetry but with line nodes at k z w ±7r/2c. We study 
this new state within both BCS and Eliashberg theories using a realistic band structure and find 
that it is more stable than the d x 2_ y 2 (Big) state over a wide range of parameters. Thus, models of 
spin-fluctuation-mediated superconductivity must consider both possibilities on an equal footing. 

PACS numbers: 74.20.-z, 74.72.-h, 74.25.Dw, 74.20.Fg 



One of the fundamental questions concerning the 
cuprate superconductors is the symmetry of the super- 
conducting gap function Ak- Experiments which probe 
the relative sign of the Josephson currents in the a and b 
directions in YBa2Cu307 (YBCO) (D are commonly con- 
sidered proof of d x 2_ y 2 pairing. However, recent theories 
which account for the CuO chains in YBCO have chal- 
lenged this interpretation 0], and other measurements 
exist 0] which are difficult to reconcile with this pairing 
symmetry. Thus, current experiments do not unambigu- 
ously establish d x 2_ y 2 pairing for all the cuprates under 
all conditions. 

The experiments of Ref. j[) are consistent with theo- 
ries of high-temperature superconductivity based on an- 
tiferromagnetic spin fluctuations Since the Cu02 
planes in the cuprates are generally considered the most 
important structural elements for the superconductivity, 
most work on this subject has focussed on a single, two- 
dimensional (2D) band, representing the physics of an 
isolated Cu02 plane. Within this framework, it is now 
well established that antiferromagnetic spin fluctuations 
yield singlet pairing in a d a .2_ J( 2-like orbital state 
This pairing symmetry can be understood by noting that 
antiferromagnetic spin fluctuations produce a repulsive 
interaction which inhibits on-site pairing but favors pair- 
ing on nearest-neighbor sites M. 

In materials like YBCO, though, there are two Cu02 
planes per unit cell. One therefore expects that two 2D 
bands will be involved in the physical description, both 
of which may contribute to the superconductivity. Re- 
cent calculations within such bilayer systems based on 
either weakly |p|-|io| or strongly 11 lj] correlated elec- 
trons concluded that an alternative to the d x 2_ y 2 (Bi g ) 
state may exist. In this state, there are two gap functions 
(one for each band) which may be strongly anisotropic in 
wave vector space but which do not have nodes; however, 
the two gap functions do have opposite signs. Hence, 
this state has been termed the s ± state [|l]|. Strong- 
coupling Eliashberg calculations have demonstrated that 
the state is actually more stable than the d x 2_ y 2 state 



whenever the inter-band antiferromagnetic correlations 
are stronger than the intra-band ones 0. This pairing 
symmetry may also be consistent with all of the ex- 
perimental determinations of the phase of Ak . 

The s ± state is restricted to the bilayer cuprates such 
as YBCO, and thus does not address the superconductiv- 
ity in single-layer cuprates like Lai_ x Sr x Cu04 (LSCO). 
Moreover, both and d x 2_ y 2 states are based on a 2D 
picture that neglects the third dimension. In particular, 
increasing the bilayer separation in a three-dimensional 
(3D) layered system causes the two-band structure to 
cross-over into a single-band one. Thus, there should be 
a one-band analog of the state in 3D which would 
compete with the d x 2_ y 2 state. 

We examine this possibility within a general model for 
antiferromagnetic spin fluctuations in a single 3D band. 
By solving both the weak-coupling BCS and strong- 
coupling Eliashberg gap equations at T c , we are able to 
compute the phase diagram for the different pairing sym- 
metries. We find that the one-band analog of the state 
exists and is stable over a large range of antiferromag- 
netic correlation lengths which parameterize our model. 
We emphasize that, even if fully oxygenated YBCO has 
d x 2_ y 2 pairing, one may be able to observe a transition to 
this s state by changing the correlation lengths through 
changes in oxygen content, pressure, and chemical sub- 
stitution. The possible experimental relevance of this 
state therefore demands that it be considered as a seri- 
ous alternative to the conventional d x 2_ y 2 state in future 
calculations. 

To be as general as possible, we take the antiferromag- 
netic spin fluctuation interaction to have the form 



V(q,u) = 



V 



1 + &(q|| - Q||) 2 + e c (q±- - Q±) 2 - iw/wsF ' 

(1) 



where Vq is the electron-spin fluctuation coupling 
strength; (QiuQj.) = (V/a, 7r/a, tt/c) and u>sf = 7.7 
meV are the spin fluctuation wave vector and charac- 
teristic frequency; and and £ c are the intra- and 
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inter-planar antiferromagnetic correlation lengths, re- 
spectively. [Throughout this paper, a (c) is the unit cell 
dimension normal to (along) the z direction.] When £ c — > 
0, Eq. ([l]) reduces to the 2D spin fluctuation model of Mil- 
lis, Monien, and Pines |l4p|]. £ c is introduced into this 
model in a straightforward way to account for the exper- 
imental fact that inter-planar antiferromagnetic correla- 
tions exist in YBCO Q. 

This interaction [Eq. (|l|)] is assumed to exist between 
quasiparticles which have a 3D dispersion relation e^. To 
further relate our calculations to the cuprates, we take 
ek from the single-band model developed by Andersen et 
al. for YBCO M: 



et = — 2t i (cos k x a + cos k y a) — At 2 cos k x a cos k y a 

— -^-(cos k x a — cos k y a) 2 cos k z c — /i. (2) 

Here, h = 0.25 eV, t 2 = -0.1 eV, and t 3 = 0.05 eV are the 
matrix elements for intra-planar nearest-neighbor, intra- 
planar next-nearest-neighbor, and inter-planar hopping, 
respectively; \i = -0.4 eV is the chemical potential. Note 
the surprising dependence of the inter-planar hopping 
on the intra-planar wave vector. This behavior results 
from the reduction of the full local-density approxima- 
tion band structure to a single-band, tight-binding dis- 
persion; see Ref. for details. Equation (||) reproduces 
the general features of the band structure of YBCO re- 
vealed by angle-resolved photoemission experiments fTijl 
and yields the Fermi surface shown in Fig. [l]. 

With the pairing interaction and the band structure 
specified, we can proceed to study the superconductivity 
within either BCS theory or its strong-coupling counter- 
part, Eliashberg theory pq| . For simplicity, let us first 
consider the BCS equation at T c , 



ApairAk — 



1 



Etanh(. 
Vic-k 

k' 



e k ,/2k B T) 
2e k < k ' 



(3) 



where Ak is the gap function, Vq = V(q, 0) is the pair- 
ing interaction [Eq. (Q)], ek is the electronic dispersion 
[Eq. (|J)], N is the number of sites in the lattice, T 
is the temperature, and A pa ir is the pairing eigenvalue. 
This equation is written for singlet pairing, for which 
A_k = Ak; for triplet pairing, A_k = — Ak, and the 
right-hand side of Eq. (|) is multiplied by -1/3. The 
Eliashberg equations are similar in spirit although for- 
mally more complicated; the interested reader is referred 
to the literature [Q. 

Equation (||) is an eigenvalue equation, allowing the 
eigenvectors Ak to be classified into different representa- 
tions of the crystal point group by their transformation 
properties under the elements of that group [H],^(J . The 
"pairing symmetry" of Ak refers to this representation. 
For example, in tetragonal systems (point group D^h), 
u d x 2_ y 2 pairing" corresponds to the B\ g representation 



FIG. 1. Gap functions Ak for a three-dimensional 
antiferromagnetic spin-fluctuation-mediated superconductor 
[Eq. ([j])] which transform under the (a) B\ g (d x 2_ y 2) and (b) 
Ai g (s ± ) representations of the tetragonal point group. The 
Fermi surface in the first Brillouin zone [Eq. (0)] is shown 
colored by the value of Ak on the Fermi surface, with blue 
denoting positive values and red denoting negative ones. Ob- 
serve the different locations of the gap nodes (white stripes) 
around the Fermi surface. 



and "s^ pairing" to the A\ g representation. To deter- 
mine the most favored pairing symmetry, one can either 
find the representation yielding the largest A pa ir at fixed 
T or the largest T c (defined by A pair (T = T c ) = 1). We 
employ both methods to construct the phase diagram for 
our model and find consistent results. 

We solve Eq. (||) and its strong-coupling counter- 
part using a fast Fourier transform technique with the 
wave vector integrations carried out over the entire Bril- 
louin zone The extremal eigenvalues for each one- 
dimensional representation of the point group are ob- 
tained by an iterative procedure with an eigenvalue shift 
p2fl . The BCS calculations are carried out on a 32 x 32 
x 32 discretization of the Brillouin zone, and the strong- 
coupling calculations are performed on a 64 x 64 x 8 
mesh with 256 Matsubara frequencies, corresponding to 
a frequency cut-off on the order of three times the band 
width. Except as noted below, our results are insensitive 
to these choices. 

In discussing our results, we begin with an exami- 
nation of the gap function symmetries associated with 
the dominant pairing eigenvalues. We have calculated 
Apair for all the one-dimensional representations of the 
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tetragonal point group over a large range of correlation 
lengths (see below), and we find only two representa- 
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tions which lead to a superconducting instability: 
and Ai g . The B\ g solution is shown in Fig. |](a) an d 
corresponds to the d x 2_ y 2-\ike solution encountered in 
previous work on a single band in 2D [Q-||. This state 
possesses nodal planes defined by the directions [110] and 
[001] and by [110] and [001], and can be thought of as an 
intra-planar pairing of quasiparticles on nearest-neighbor 
sites 0. On the other hand, the A\ g solution is new. 
From Fig. |(b), we see that the A± g solution, while for- 
mally having "s-wave" (i.e., A± g ) symmetry, possesses 
nodal lines at k z as ±7r/2c j23|. The gap function Ak 
thus transforms like Aj. ~ cos k z c and corresponds to the 
pairing of quasiparticles between different planes. This 
solution is the one-band analog of the s ± state discussed 
within the context of two-band, 2D models |8|-|l^|. We 
can see this analogy by mentally translating every other 
plane in this monolayer model to produce a bilayer; under 
this transformation, the regions of Ak with the same sign 
in Fig. 0(b) would be mapped onto one of the resulting 
two bands, and the nodes would disappear in favor of the 
band gap, maximizing the superconducting condensation 
energy. 

The one-band s state exists as a possible supercon- 
ducting solution, but is it stable? When £ c — > 0, Eq. (||) 
indicates that Ak can only depend on k x and k y ; we 
therefore expect that the only stable solution will be the 
d x 2 _y2-\ihe (Big) one. On the other hand, when £ t — * 0, 
Eq. (H) must yield a Ak which depends only on k z . More- 
over, we require a sign change in the Ak connected by 
wave vectors where Vq is maximal (at q = (0, 0, tt/c)) in 
order to cancel the negative sign on the right-hand side 
of Eq. (S). These demands are met by the s ± state of 
Fig. 0(b), and so this solution should be stable in this 
limit. For intermediate values of the correlation lengths, 
we should see a cross-over between these states [M . 

We have tested these expectations numerically and 
present our results in Figs. ||and||. Consider first Fig. ||, 
which shows A pa i r as a function of £ c for £ a f, = 2.3a 
computed within strong-coupling Eliashberg theory. The 
largest pairing eigenvalues, and hence the most favored 
superconducting solutions, belong to the B± g (d x 2_ y v) or 
A\ g [s ) representations. As £ c increases, we see that 
the most stable solution crosses over from the B\ g to the 
A\ g representation, as we expect from the preceding dis- 
cussion. 

To give these results a more physical interpretation, 
we show in Fig. || the critical temperature T c of the most 
stable pairing symmetry as a function of the intra- and 
inter-planar correlation lengths determined from a BCS 
calculation p5 26 1. The A± g state is stable over a large 



region of parameter space, with the phase boundary given 
roughly by the line £ c /c = 2 £ a (,/a. The limiting cases 
discussed above are confirmed as well: the B\ g (d x 2_ y i) 
solution dominates as £ c — ► 0, and the reverse is true for 
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FIG. 2. Pairing eigenvalues A pa ir belonging to the B\ g 
(d x 2_ y 2, squares) and Ai g (s , triangles) representations as 
a function of the normalized inter-planar antiferromagnetic 
correlation length £ c /c; the solid lines are to guide the eye. 
These results are computed within strong-coupling Eliashberg 
theory for an antiferromagnetic spin-fluctuation-mediated in- 
teraction [Eq. (|l|)] with an intra-planar correlation length £ a & 
= 2.3a at 90 K. Note the crossover between Bi g and Ai g pair- 
ing around £ c = 3.5c, which is consistent with the BCS result 
of Fig. | 



the Ai g (s^) state as £ (, — > 0. In addition, the critical 
temperatures are largest along the £ c = and £ a b = 
lines and increase as the antiferromagnetic correlation 
lengths increase. 

Moving towards the phase boundary in Fig. ||, T c de- 
creases but does not vanish, at least for large £ a h and £ c . 
At small correlation lengths, however, T c is small near 
the phase boundary and very sensitive to the finite-size 
effects induced by our discretization of the Brillouin zone. 
Hence, we cannot exclude the possibility of a finger of 
normal phase extending between the two superconduct- 
ing phases at small £ a b and £ c . With this caveat, we 
believe the Fig. || represents the true phase diagram of 
our model. 

It is important to note that the existence of the one- 
band s (Ai g ) state depends only on the 3D interaction 
and not the detailed band structure. Specifically, recom- 
puting the phase diagram with t-j = in Eq. (^) gives 
nearly identical results. The s state may therefore be 
compatible with theories of the cuprate superconductors 
in which single-particle inter-layer hopping is suppressed 
but Cooper pair tunneling is not |^7j . Of course, the sta- 
bility of the state in any concrete system, as opposed 
to its existence, depends on microscopic details such as 
the band structure and must therefore be evaluated on a 
case-by-case basis. 

To conclude, we have shown that single-band models 
of antiferromagnetic spin fluctuation-mediated supercon- 
ductivity in 3D do not automatically give rise to d x 2_ y 2 
pairing. Rather, a phase which possesses formal "s-wave" 
(A\g) symmetry with nodal lines at k z w ±7r/2c exists 
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FIG. 3. BCS critical temperature T c of an antiferromag- 
netic spin-fluctuation-mediated superconductor [Eqs. 
as a function of the normalized inter- (£ c /c) and intra-planar 
(£ab/a) correlation lengths. Light shading represents T c for 
Bi g {d x 2_ y 2) pairing while dark shading represents Ai g (s ) 
pairing. 



and is stable over a wide range of parameters. Within 
the context of the cuprate superconductors, these results 
imply that this unusual s state, which is the one-band 
analog of the state of Ref. Q , must be considered on 
equal footing with the conventional d x 2_ y 2 state. This 
point is reinforced by the fact that the s state is fa- 
vored when the inter-planar antiferromagnetic correla- 
tions are stronger than the intra-planar ones, a situa- 
tion which may obtain in the cuprates. Finally, the re- 
sults of this paper suggest that, even in a superconductor 
with pure d x 2_ y 2 pairing, substitution, doping, or pres- 
sure may lead to a change in the correlations lengths 
and thus to a transition between different gap function 
symmetries. Such transitions may already have been ob- 
served in angle-resolved photoemission experiments p8[ . 
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